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ABSTRACT 

A  computer  mathematical  model  that  will  simulate  the  six  degree 
of  freedom  motion  of  a  submarine  has  been  developed.  Hydrodynamic 
coefficients  obtained  from  the  testing  of  physical  models  enable  the 
user  to  accurately  simulate  specific  submarine  designs.  The  simulation 
model  can  be  used  as  a  design  tool  to  study  and  predict  submarine 
dynamic  response.  Complete  three  dimens  tonal  motion  1i  allowed  for 
the  submarine  while  constraints  In  powering,  rudder  deflection,  dive 
plane  ingle,  and  response  time  can  be  selected  by  the  use'.  Outputs 
Include  a  chronological  history  of  all  linear  and  angular  velocities 
and  accelerations,  the  Instantaneous  angle  of  deflection  of  <?ach 
control  surface,  and  the  trajectory  of  the  submarine  center  of 
gravl ty.  _ 
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ABSTRACT 

A  computer  mathematical  model  that  will  simulate  the  six  degree 
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coefficients  obtained  from  tne  testing  of  physical  models  enable  the 
user  to  accurately  simulate  specific  submarine  designs.  The  simulation 
model  can  be  used  as  a  design  tool  to  study  and  predict  submarine 
dynamic  response.  Complete  three  dimensional  motion  Is  allowed  for 
the  submarine  while  constraints  In  powering,  rudder  deflection,  dive 
plane  angle,  and  response  time  can  tie  selected  by  the  user.  Outputs 
Include  a  chronological  history  of  all  linear  and  angular  velocities 
and  accelerations,  tne  instantaneous  angle  of  deflection  of  each 
control  surface,  and  the  trajectory  of  the  submarine  center  of 
grsvl ty . 
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CHAPTER  I  -  INTRODUCTION 


I.}  background 

During  the  several  \ ears  spent  in  designing  a  submarine,  a 
great  deal  of  effort  Is  expended  trying  to  Imp-ove  the  design  so  that 
the  best  possible  vessel  Is  eventually  constructed.  During  the  de¬ 
sign  phases  tne  designer  will  frequently  examine  problems  that  were 
found  in  older  submarines  and  try  to  determine  the  cause  of  each 
problem  so  that  they  can  be  eliminated  In  the  new  design.  At  several 
stages  during  the  design,  computer  models  are  used  to  assist  In  such 
things  as  structural  design,  equilibria  calculations,  and  internal 
arrangement.  As  the  design  progresses,  physical  models  of  the 
proposed  vessel  are  built  and  tested  In  a  towing  tank  to  measure  the 
nydrodynamlc  coefficients  finally,  when  It  Is  possible  to  deal  with 
the  vessel  as  a  whole,  a  dynamic  analysis  Is  conducted  to  gauge  how 
the  submarine  will  perform  in  an  underwater  environment  when  ill  six 
degrees  of  freedom  are  available. 

Wiile  it  is  possible  to  gather  dynamic  information,  such  as  the 
hydrodynamic  coefficients ,  from  physical  model  tests,  the  models  are 
necessarily  too  restricted  in  their  motion  to  permit  a  complete 
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analysis.  It  Is  necessary  to  construct  a  mathematical  model  of  the 
submarine  to  simulate  the  submarine's  motion  and  thereby  obtain 
•nou^n  data  for  a  complete  analysis.  The  simulation  model  will  use 
the  hydrodynamic  coefficients  obtained  from  the  phyilcal  model  as  a 
basis  for  the  simulation*  A  properly  constructed  model  will  permit 
the  designer  to  simulate  any  conceivable  maneuver  and  gauge  the 
submarine's  response.  The  simulation  model  then  becomes  a  tool  to 
assist  In  improving  the  design  and  achieving  the  best  possible 
vessel.  After  the  submarine  Is  designed,  the  model  can  continue  to 
be  useful  by  assisting  In  the  evaluation  of  operating  procedures. 

The  model  can  be  placed  In  any  maneuvering  situation  without  hazard 
to  crew  or  vessel.  This  Is  especially  useful  In  evaluating  casualty 
situations.  The  model  can  also  be  used  to  estimate  the  effect  of 
proposed  design  changes  to  the  submarine.  For  Instance,  the  model 
can  snow  the  effect  of  changing  the  maximum  deflection  of  a  control 
su-face  cr  its  rate  of  operation. 

All  of  the  few  six  degree  of  freedom  submarine  simulation  models 
In  existence  are  too  expensive  for  dally  use  in  the  design  office. 
Gunenl  1y,  the  models  are  complicated  to  use  and  sometimes  difficult 
to  obtain.  These  problems  nave  created  the  need  for  a  program  that 
Is  both  Inexpensive  enougn  to  permit  dally  use  In  the  design  office 
and  simple  enough  so  that  anyone  can  use  it.  L.is  need  has  formed 
the  motivation  for  this  thesis. 

1.2  MQOEL  DEVELOPMENT 

A  mathematical  model  must  use  Newton's  law  of  Motion  and  the 
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differential  equations  resulting  from  a  dynamic  analysis  of  the 
submarine.  A  submarine  has  movable  appendages,  so  these  must  also 
be  accounted  for  In  the  model.  Newton’s  Law  of  Motion  can  be 
expressed  as: 

(1)  Force  *  ^  (momenta* ) 

(2)  Moment  ■  ^  (angular  momentum) 

If  a  submarine  with  a  coordinate  system  such  as  shown  In  Figure  1  Is 
used,  then  equation  (1)  can  be  applied  along  each  of  the  three  axes. 
Similarly,  equation  (2)  can  be  applied  around  each  axis.  This  gives 
a  total  of  six  equations  to  represent  the  six  degrees  of  freedom  of 
the  submarine.  These  six  equations  of  motion  are  well-known  [1]  so 
they  are  not  developed  In  detail  here.  They  are,  however.  Included 
for  reference  In  Appendix  A. 

fchlle  the  equations  of  motion  for  the  submarine  form  the  heart 
of  the  model,  they  are  Insufficient  by  themselves.  On  an  actual 
submarine,  the  officer  of  the  deck  orders  the  rudder  or  dive  planes 
moved  In  order  to  maneuver  the  ship.  The  routines  will  normally 
sense  where  the  ship  currently  is  and  where  it  Is  supposed  to  be  and 
then  move  the  appendages  in  the  appropriate  direction  just  as  the 
deck  officer  would  do.  The  six  equations  of  motion  together  with 
the  appendage  control  subroutines  constitute  the  vital  components  cf 
the  simulation  model.  A  main  program  Is  necessary  to  coordinate  the 
Input,  output,  and  to  control  the  action  of  the  subroutines. 

The  components  cf  the  model  are  discussed  In  detail  in  Chapter 
II  of  this  thesis.  The  model  Is  written  as  a  computer  program  whose 
features  are  discussed  In  Chapter  III.  The  final  chapter  will 


discuss  the  model  tests  and  their  results.  A  listing  of  the 
program  will  be  provided  In  Appendix  B. 


12 


CHAPTER  II  -  THE  SIMULATION  MODEL 

II. 1  GENERAL  CAPABILITIES 

The  simulation  model  dev*  1  ope o  In  this  thesis  will  provide 
trajectory  Information  for  routine  submarine  maneuvers  such  as 
turning  or  changing  depth.  The  trajectory,  which  Is  referenced  to 
a  fixed  coordinate  system.  Is  computed  as  a  function  of  time.  As 
the  trajectory  of  the  model  Is  developed,  the  velocities  and 
accelerations  are  calculated  and  stored  for  output  to  the  user. 

Both  angular  and  linear  velocities  and  accelerations  are  provided. 

The  angle  of  deflection  of  each  control  surface  appendage  is  computed 
for  each  step  of  the  entire  trajectory  or  maneuver. 

The  rate  of  deflection  of  each  control  surface  and  Its  angle 
of  maximum  deflection  are  specified  by  the  user.  Any  of  the  control 
surfaces  can  be  "jaamed*  by  specifying  a  particular  angle  of 
deflection.  The  model  can  be  Initially  placed  on  any  course  and 
depth  and  then  ordered  to  come  to  any  new  course  and  depth.  The 
user  can  select  the  Initial  speed  and  the  appropriate  thrust 
coefficients  to  cause  a  cnange  In  speed. 


II. 2  COORDINATE  SYSTEMS 
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The  coordinate  system  plan  used  In  the  simulation  model  Is 
based  on  [1].  The  model  uses  two  orthogonal  coordinate  systems;  one 
remaining  fixed  at  the  water  surface  while  the  otheh  travels  with  the 
submarine  to  act  as  a  local  reference  system.  The  fixed  system, 
designated  xQ,  yQ,  z0,  Is  related  to  the  moving  system,  designated 
x*  y.  x,  through  the  angles  t>  ,  *  as  Illustrated  In  Figure  2. 
Quantities  measured  In  the  moving  system  can  be  referenced  to  the 
fixed  system  by  using  the  transformation  [2]  given  In  Appendix  A. 

The  origin  of  the  moving  system  Is  at  the  center  of  gravity  of 
the  submarine.  This  has  the  advantage  that  only  the  principal 
moments  of  Inertia,  Ix,  l^,  I2,  are  non-zero  (l.e.:  Ixy“IyZ*Ix2"  °)- 
It  has  the  additional  advantage  that  the  equations  of  motion  In  [1] 
also  use  this  reference  point  for  the  x,  y,  *  system.  It  would  have 
been  possible  to  use  the  centerline  of  the  submarine  for  the  origin. 
This  would  have  the  advantage  of  making  better  use  of  the  symmetry 
of  the  vessel,  howevur  It  would  nave  the  disadvantage  of  having  to 
correct  both  the  equations  of  motion  and  some  of  the  hydrodynamic 
coefficients  for  the  new  origin. 

Appendage  movements  are  measured  with  respect  to  the  moving 
coordinate  system.  All  velocities  and  accelerations  are  measured 
along  the  axes  of  the  movl  .g  system.  This  Is  also  shown  In  Figure  2. 

II. 3  EQUATIONS  QF  NOTION 

The  general  nature  of  the  six  degree  of  freedom  model  requires 
the  use  of  the  six  equations  in  their  non-linear  form.  A  large 
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nustoer  of  hydrodynamic  coefficients  are  necessary  for  this  model. 

It  would  be  desirable  to  derive  each  coefficient  from  existing  hydro¬ 
dynamic  theory.  Some  coefficients  have  been  obtelned  for  todies  of 
revolution  or  other  mathematically  amenable  shapes however,  when 
appendages  such  as  control  surfaces,  falrwaters,  and  propellers  are 
taken  into  account,  the  accuracy  of  the  theoretical  values  Is  brought 
Into  question.  For  this  reason  It  Is  standard  practice  to  obtain  the 
ntmrical  value  of  the  coefficients  by  experimental  means.  This 
usually  entails  the  towing  and  measurement  of  physical  models. 

Once  the  hydrodynamic  coefficients  are  known,  the  equations 
can  be  solved  for  the  accelerations.  The  six  equations  of  motion 
must  be  written  such  that  the  highest  order  derivatives,  namely 
u,  v,  w,  p,  q,  r  and  their  coefficients ,  appear  on  the  left  hand 
side  of  the  equation.  This  gives  each  equation  a  form  like: 


*  *  ‘1,4.2  “  *  *1,4*3  0  *  *1.4*4  "  *  *1.4*5  r 


*1.4  “  *  *1.4*! 

■  fj  (u,  v,  w,  p,  q,  r,  o,  ,  jr,  j^.,  5^»  *  1,  m,  .  .  .  ) 

The  left  hand  side  Is  placed  In  a  matrix  format.  The  six  equations 
are  then  represented  as: 
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The  method  of  solution  Is  an  Iterative  technique  In  which  Initial 
values  are  used  In  the  right  hand  side  functions  and  the  six  acceler¬ 
ations  are  solved  for  on  the  left  side.  The  accelerations  are  then 
used  tc  update  the  right  side  functions.  A  small  time  Increment  Is 
made  and  tne  accelerations  are  again  solved  by  using  the  latest 
value  for  the  right  side  functions. 

With  each  Iteration  the  accelerations  are  used  In  a  Taylor 
series  expansion  to  calculate  the  velocities.  The  calculations 
have  the  form: 


( t>At)* 


(t)  ♦ 


*Ut) 


The  pitch,  roll,  and  yaw  ar  jles  are  calculated  In  a  similar  manner 
using  the  angular  velocity  and  acceleration. 


(t  ♦  st)  • 


(t)  ♦ 


(St)  * 


■Of) 


In  order  to  plot  the  trajectory  of  the  submarine,  a  coordinate 

system  transformation  Is  performed  to  obtain  the  linear  velocities 
•  •  • 

x  .  yQ ,  zQ  Hi  the  fixed  coordinate  system.  The  velocities  are  used 
In  a  Taylor  series  expansion  to  obtain  the  position  of  the  center  of 
gravity  of  the  submarine  In  the  xQ,  yQt  zQ  system. 
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When  the  position  of  the  submarine  Is  calculated  the  time  Is  advanced 

one  Increment  and  another  Iteration  begins. 

The  Iteration  loop  will  continue  to  operate  until  so ue  test 

criteria  are  met.  For  Instance.  If  the  model  Is  performing  a  dive. 

then  the  model  tests  the  value  of  z„  to  see  If  the  model  Is  at  the 

o 

appropriate  depth.  A  test  Is  also  made  of  the  pitch  angle  to  see  If 
It  Is  within  son*  specified  range  of  zero,  lastly  a  check  Is  made  to 
ensure  that  all  of  thi  dive  planes  are  at  zero  deflection. 

II. 4  CONTROL  SURFACES 

The  motion  of  the  model  Is  controlled  by  the  movement  of  the 
control  surfaces.  The  control  surfaces  Include  the  rudder,  the 
stem  planes,  and  the  sail  planes.  For  analysis  of  submarine  designs 
the  movement  of  the  control  surfaces  can  be  governed  by  an  automatic 
control  system.  In  the  model  each  control  surface  Is  provided  with 
Its  own  automatic  control  system.  The  principle  of  operation  for 
all  of  the  automatic  control  systems  Is  the  same.  In  each  case  the 
deflection  of  the  control  surface  Is  made  proportional  to  an  error 
signal  and  the  rate  of  change  of  the  error.  The  sail  planes,  for 
Instance,  control  the  depth  of  the  model.  The  calculated  deflection 
of  the  planes  Is  proportional  ,,o  the  depth  error  and  the  rate  of 
change  of  depth. 
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*c  •  (present  depth  -  ordered  depth) 

♦  *2  (rate  of  change  of  depth) 

A  large  error  will  Initially  cause  a  large  deflection,  but  as  the 
rate  of  change  Increases,  the  deflection  will  decrease  until  some 
moderate  rate  Is  achieved.  Since  movement  of  the  sailplanes  does 
not  have  any  significant  effect  on  any  motion  other  than  depth. 

Its  control  system  uses  only  the  variables  associated  with  depth. 

The  stem  planes  control  pitch  angle  as  well  as  depth.  The 
control  systom  accounts  *or  this  by  using  the  pitch  angle,  <7,  and  the 
rate  of  change  of  4. 

Sc  “  *3  (?rw*nt  depth  -  ordered  depth) 

♦  (rate  of  change  of  depth) 

♦  (pitch  angle) 

♦  Kg  (rate  of  change  of  pitch  angle) 

The  K's  In  the  control  systems  are  proportionality  constants 
known  as  gain.  Tie  value  of  the  gain  determines  the  sensitivity  and 
response  of  the  control  system.  Since  the  error  signal  continually 
varies  during  a  maneuver.  It  Is  best  to  keep  the  gain  low.  It  Is 
desirable  to  have  just  enough  gain  on  the  error  signal  to  get  the 
control  system  moving  and  keep  it  moving  In  the  right  direction. 

The  gain  on  tne  rate  signals  should  be  much  stronger  to  dampen  out 
the  oscillations  caused  oy  response  to  the  error  signal. 

The  control  surfaces  move  at  a  rate  specified  by  the 
user.  Whenever  the  calculated  deflection  differs  fro v*  the 
present  deflection,  the  control  surface  will  move  toward 
the  calculated  value  at  the  specified  rate. 
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CHAPTER  THREE  -  PROGRAM  USER’S  GUIDE 

111. 1  SENERAL  FEATURES 

The  simulation  nodal  consists  of  a  main  program  and  four 
subroutines.  The  four  subroutines  are  named  FUNC.  RUDDER,  DEPTH, 
and  STERN.  The  main  program  reads  the  Input  data,  changes  units, 
Initializes  values,  and  prints  the  output.  The  main  program  performs 
the  Taylor  series  expansions,  the  coordinate  system  transforms,  and 
provides  the  logical  statements  for  calling  the  subroutines.  On  the 
first  Iteration,  when  time  equals  zero,  no  subroutines  are  called; 
the  output  from  this  first  Iteration  is  then  a  statement  of  the 
Initial  conditions  of  the  problem.  Each  succeeding  Iteration  will 
call  subroutine  FUNC  to  calculate  the  new  position  of  the  model. 

Calls  to  the  control  surface  subroutines  are  made  on  the  basis  of 
need,  with  no  subroutine  being  called  more  than  once  In  the  Iteration. 
The  decision  on  whether  or  not  to  call  RUDOER,  DEPTH,  or  STERN  Is 
contained  In  a  series  of  logical  IF  statements  In  the  main  program. 

1 1 1. 2  SUBROUTINE  FUNC 


Subroutine  FUNC  contains  the  six  equations  of  motion.  The 
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Input  parameters  to  FUNC  consist  of  the  hydrodynamic  coefficients, 
the  (i  ,del  **eloc1t1es,  the  orientation  «r*  space,  all  of  the  ship's 
characteristics  such  as  length  and  moments  of  Inertia,  the  deflection 
of  each  control  surface,  and  the  propulsive  coefficients.  The  sub¬ 
routine  returns  the  value  of  the  six  accelerations:  UOT,  VOT,  WOT, 
POT,  QOT,  ROT.  The  Input  and  output  for  FUNC  are  all  contained  In 
CQfflON  /FOUR/  and  COWON  /FIVE/.  The  solution  to  the  matrix  equation 
In  subroutine  FUNC  Is  made  possible  by  a  call  to  the  library  function 
LEQT1F  [4].  LEQT1F  performs  a  Gaussian  reduction  for  the  matrix 
equation.  Any  similar  Gaussian  reduction  could  be  substituted  If 
LEQT1F  Is  not  available.  An  explanation  of  the  parameters  used  In 
the  call  to  LEQT1F  Is  found  In  Appendix  B. 

111. 3  SUBROUTINE  RUOOER 

Subroutine  RUOOER  control-  ail  horizontal  motion  for  the 
model.  The  input  parameters  Include:  K9,  K1Q,  T5,  T6,  Til,  T12, 

T17,  T18,  TLA6R,  RRATE,  RUOAKT,  COURSE.  OaR,  DELT,  R.  PSI.  The  sub¬ 
routine  :etun.s  a  new  value  for  DEIR,  the  rudder  deflection,  based  on 
a  calculation  using  its  present  and  ordered  headings.  The  Inputs  are 
all  contained  In  COWON  /THREE/  and  COr’ON  /FIVE/ .  Subroutine  RUOOER 
Is  called  whenever  the  model  Is  not  on  the  desired  course  or  when  the 
rudder  Is  deflected. 

111. 4  SUBROUTINE  DEPTH 

Subroutine  DEPTH  controls  the  forward  set  of  dive  planes.  It 
can  control  either  bow  planes  or  sail  planes  depending  on  what  the 
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submarine  Is  fitted  with.  This  subroutine  Is  sensitive  to  the  depth 
error  end  the  rate  of  change  of  depth.  Subroutine  DEPTH  will  move 
the  forward  planes  In  the  direction  necessary  to  bring  the  depth 
error  to  zero.  As  Input  parameters,  the  subroutine  uses:  K1 ,  K2 
Tl,  T2 ,  T13,  T14,  T15,  T16,  TLAG3,  DELT,  DIFF,  AQIFk,  2DT,  ATHETA, 
MAXANG,  DCR1T.  The  subroutine  returns  a  new  value  for  DELB,  the 
bow  plane  deflection,  after  each  call. 

If  MAXANG,  the  maximum  dive/ ascent  angle,  has  been  exceeded, 
then  DEPTH  will  return  a  diagnostic  write  statement.  The  user  may 
speclty  MAXANG,  the  maxlimja  pitch  angle  for  the  model.  Whenever 
MAXANG  Is  exceeded  the  diving  planes  are  moved  to  reduce  the  pitch 
angle  and  a  diagnostic  signal  is  generated  from  within  subroutine 
DEPTH.  The  diagnostic  will  say  "Maxlmim  dive/ascent  angle  exceeded 

at  time  _ .  Standard  fixup  taken.-  Since  the  dive  planes  only 

begin  to  react  when  MAXANG  is  exceeded,  the  submarine  will  oversnoot 
the  angle  before  a  reduction  in  the  angle  actually  occurs.  The 
diving  planes  will  continue  to  operate  until  the  pitch  angle  Is 
less  than  MAXANG. 

DEPTH  makes  only  a  minimal  attempt  to  control  the  pitch  angle 
of  the  model;  it  only  warns  the  user  wnen  the  angle  Is  exceeded,  and 
then  moves  the  bow  planes  so  tnat  the  pitch  angle  does  not  grow 
larger.  The  principal  purpose  of  DEPTH  is  to  provide  depth  control; 
it  does  this  In  conjunction  with  subroutine  STERN. 

III. 5  SUBROUTINE  STERN 


Subroutine  STERN  controls  the  movement  of  tne  stem  planet  to 
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achieve  the  rieslnsd  depth  and  pitch  angle.  The  movement  of  the 
planes  Is  sensitive  to  the  depth  error ,  the  rate  of  change  of  depth, 
the  pitch  angle,  and  the  rate  of  change  of  pitch.  The  stem  planes 
Mill  move  to  bring  the  depth  error  to  zero  and  the  pitch  angle  to 
zero.  As  Input  parameters  the  subroutine  uses:  KS,  K6,  K7,  K8,  T3, 
T7,  T8,  T10.  TLAGS,  STERAT,  STLRMX,  THETA,  Q,  DELS.  DIFF ,  ADIFF, 

ZDT,  ATHETA,  MAXANG,  NOPICH,  DCRIT.  The  subroutine  returns  a  new 
value  for  the  stem  plane  deflection,  DELS.  The  Input  and  output 
parameters  are  all  contained  In  COMMON  /ONE/,  COMMON  /FIVE/,  and 
COMMON  /SIX/.  STERN  Is  called  to  achieve  a  depth  change  or  to  correct 
the  pitch  angle. 

III. 6  RANGE  VARIABLES 

Since  It  is  not  usually  possible  to  bring  a  computer  simulation 
model  to  a  precise  depth  or  angle.  It  Is  necessary  to  define  ranges 
aroimd  the  desired  depth  or  angle  which  will  be  acceptable  to  the 
user.  For  example.  If  the  submarine  were  to  make  a  depth  change  of 
500  feet,  the  dive  may  be  considered  complete  if  the  model  settled 
within  ten  feet  of  the  desired  depth.  The  range  is  specified  by  the 
user  to  enable  him  to  achieve  wnatever  precision  Is  desired.  Within 
this  range  the  main  program  will  not  make  a  call  to  the  control  sur¬ 
face  subroutine',  since  no  call  Is  made,  the  control  surfaces  cannot 
be  moved.  In  the  program,  the  variable  name  for  the  depth  range  is 
NQOIFF;  It  Is  usually  set  at  ♦  5  or  ♦  10  feet.  Subroutine  DEPTH  will 
be  called  whenever  the  model  is  off  of  Its  ordered  depth  by  more  than 
NOOIFF.  The  range  about  zero  pitch  angle  is  given  the  name  NOPICH; 
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It  Is  usually  set  et  ►  1  degree.  STERN  Is  called  If  the  pitch  angle 
Is  greater  than  NOPICH.  The  variable  name  for  achieving  tne  proper 
course  Is  ONCRS;  It  Is  usually  set  at  ♦  1  degree.  RUDDER  will  be 
called  whenever  the  model  Is  off  course  by  more  than  ONCRS. 

An  additional  range  variable  Is  used  during  diving  or  surfacing. 
When  the  model  moves  within  some  critical  range  of  the  desired  depth. 

It  Is  tine  to  begin  moving  the  dive  planes  to  zero  angle  and  let  the 
vessel  glide  into  the  desired  depth.  This  maneuver  prevents 
oscillation  of  the  dive  planes  as  the  error  signal  and  the  rate 
signal  both  become  snail.  This  also  ensures  that  the  planes  are  at 
or  near  zero  angle  by  the  tine  the  desired  depth  Is  reached.  The 

m 

name  of  this  variable  is  DCRIT;  it  is  usually  set  at  50,  75  or  100 
feet  depending  on  the  speed  of  the  submarine.  STERN  will  be  called 
whenever  the  model  is  off  its  ordered  depth  by  more  than  DCRIT. 

III. 7  PR06RAH  PARAMETERS 

The  Iterative  nature  of  tne  program  relies  on  a  time  Increment 
being  made  after  each  step.  The  size  of  the  time  increment  is 
optional;  the  smaller  the  increment  the  more  accurate  the  calculations. 
In  selecting  the  size  of  the  time  Increment,  one  must  bear  <n  mind 
the  length  of  time  required  to  complete  the  maneuver  and  the  storage 
capacity  of  the  program,  jue  to  the  quantity  of  information  that  is 
calculated  by  the  program,  it  Is  practical  to  print  only  half  of  the 
data  at  one  time.  The  other  half  of  the  data  is  stored  in  the 

array  STOW.  This  data  Is  then  printed  wnen  the  maneuver  Is  complete. 
STOW  nas  the  capacity  to  hold  the  data  resulting  from  600  iterations. 
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A  time  check  Is  Incorporated  In  the  main  program  to  enable  the 
user  to  limit  the  nunber  of  Iterations.  Achieving  a  nunber  of 
iterations  equal  to  the  variable  ICNT  will  cause  the  pmiram  to  stop 
the  present  maneuver,  print  out  all  data,  and  see  If  another  maneuver 
is  desired.  ICNT  should  not  be  made  larger  than  60 0  so  that  the 
available  storage  space  is  not  exceeded. 

Each  control  surface  subroutine  has  a  time  lag  scheme  which 
senses  the  initial  conaand  to  the  control  surface  and  prevents 
immediate  action.  A  time  lag  occurs  each  time  the  direction  of 
movement  is  changed.  A  different  time  lag  may  be  specified  for 
each  control  surface. 

The  variable  INOEX  is  used  to  allow  the  user  to  perform  more 
than  one  maneuver  with  a  given  submarine  and  then  shift  submarines 
to  perform  more  maneuvers.  If  INOEX  is  less  than  or  equal  to  zero 
then  the  same  submarine  coefficients  will  be  used  for  each  set  of 
Irltlal  conditions.  If  INOEX  is  greater  than  zero  the  program  will 
reau  new  coefficients  as  well  as  a  new  set  of  initial  conditions. 


CHAPTER  FOUR  -  TEST  RESULTS  AND  CONCLUSIONS 


IV. 1  FUNDAMENTAL  ACTION  TEST 

Thu  validation  of  the  simulation  model  Is  a  matter  of  Impor¬ 
tance.  The  method  used  was  an  Independent  test  of  each  component  of 
the  model,  and  then  a  series  of  tests  w1*n  the  components  of  the 
model  working  together.  The  test  maneuvers  were  selected  either 
because  the  correct  dynamic  response  was  known  or  because  the  maneuver 
was  simple  enough  so  that  the  general  nature  of  the  response  could 
be  predicted. 

The  first  portion  of  the  model  to  be  tested  was  the  subroutine 
FUNC.  FUNC  obtains  the  solution  to  the  six  equations  of  motion.  It 
was  Important  to  establish  that  these  equations  were  properly 
Installed  In  the  program.  A  test  of  subroutine  FUNC  necessitated  a 
simultaneous  test  of  the  MAIN  program  to  read  in  the  hydrodynamic 
coefficients,  set  up  the  "A"  matrix  for  FUNC,  perform  the  Taylor 
series  expansions,  and  write  out  the  results.  The  test  for  FUNC 
was  to  reduce  the  GM  of  the  vessel  to  zero  and  then  to  assime  a 
constant  angle  on  the  dive  planes.  If  the  model  was  working  cor¬ 
rectly  It  should  traverse  a  perfect  circle  in  tne  vertical  plane. 
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It  mas  known  that  the  rudder  should  not  move  and  the  model  should 
not  roll  or  yaw.  The  model  should  assume  a  circular  trajectory  with 
a  constant  angular  velocity  and  a  constant  vertical  component  of 
velocity,  w.  Figure  3  shows  the  results  of  the  maneuver.  A  circular 
trajectory  was  quickly  achieved.  The  angular  velocity  was  constant, 
w  was  constant,  and  the  position  of  the  model  In  the  fixed  coordinate 
system  confirmed  the  circular  path.  The  modal  did  not  roll  or  yaw 
and  the  rudder  did  not  move. 

IV. 2  HORIZONTAL  MOTION  TEST 

With  the  knowledge  that  the  equations  of  motion,  the  Taylor 
series  expansions,  and  the  coordinate  system  transforms  are  operating 
properly,  the  subroutine  RUDDER  was  the  next  component  to  test.  It 
was  decided  that  a  simple  left  turn  of  40  degrees  would  be  an 
appropriate  test.  A  left  turn  was  chosen  Decause  that  presents  the 
greatest  opportunity  for  error.  The  original  course  would  be  000 
degr«ts  true  and  the  new  course  would  be  320  degrees  true.  Since  the 
yaw  angle  Is  positive  when  turning  rignt,  the  model  must  cope  with  a 
negative  yaw  angle  as  well  as  a  proper  method  for  dealing  with  course 
neadlngs  given  In  true  bearing.  The  subroutines  DEPTH  ana  STERN  were 
rendered  Inoperative  to  give  the  opportunity  of  seeing  RUDDER  operate 
without  Interference.  The  model  would  be  checked  to  ensure  that  the 
proper  rudder  rate  was  used  ana  that  the  rudder  angle  did  not  exceed 
the  maxi men  angle  ordered.  The  model  would  be  expected  to  roll  and 
squat  In  the  turn.  A  change  In  depth  was  expected  since  the  dive 
planes  could  not  act.  As  the  model  approached  the  new  heading  the 
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rudder  should  bo  put  amidships  and  the  vessel  should  steady  out. 

In  Figure  4  the  relationship  between  the  rudder  and  the  heading 
Is  shown  as  a  function  of  time.  Tne  model  reacted  as  predicted 
exct  t  possibly  at  the  end  of  the  turn.  The  rudder' was  amidships  and 
the  new  heading  was  achieved  but  the  program  did  not  run  long  enough 
to  ensure  that  the  model  was  steady  on  the  course.  The  model  did  roll, 
squat,  and  change  depth  as  predicted.  The  test  for  RUDOER  was 
considered  to  be  accurate  and  sufficiently  complete  to  warrant  moving 
to  the  next  test. 

IV. 3  VERTICAL  MOTION  TEST 

Since  DEPTH  and  STERN  operate  together  to  regulate  the  depth 
and  pitch  of  the  model,  they  were  tested  together.  The  test  consisted 
of  a  simple  dive  with  a  depth  change  of  700  feet.  The  course  was  not 
changed  so  the  rudder  should  not  move.  The  model  should  not  roll  or 
yaw.  The  dive  planes  should  move  in  the  proper  direction  and  at  the 
ordered  rate.  They  should  not  deflect  more  than  the  ordered  angle. 

The  model  should  pitch  downward  and  If  the  maximum  ordered  pitch 
angle  Is  exceeded  then  the  dive  planes  should  move  to  reduce  the 
pitch  angle.  As  It  approaches  the  desired  depth,  the  model  should 
slow  Its  rata  of  descent  and  settle  within  tan  feet  of  the  depth, 
with  0  -1  degree  of  pltcn  angle. 

The  trajectory  for  this  test  Is  shown  In  Figure  5.  The  model 
acnleved  steady  state  only  six  feet  beyond  the  desired  depth;  the 
pitch  angle  was  >.03  degrees.  The  dive  planes  were  all  at  zero 
deflection.  The  model  stayed  within  ten  feet  of  the  ordered  depth 
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and  within  one  degree  of  zero  pitch  angle  for  twenty  seconds.  This 
was  done  to  ensure  that  a  steady  state  had  been  achieved. 

IV. 4  COMPLETE  HOPEI.  TEST 

As  a  final  test  the  entire  simulation  model  must  work  together. 
This  test  repeated  the  forty  degree  left  hand  turn  but  this  time  the 
dive  planes  were  allowed  to  mact  to  try  to  maintain  the  depth.  The 
program  was  kept  running  until  th3  pitch  angle  was  within  one  degree 
of  being  zero,  the  model  was  within  one  degree  of  the  proper  course 
and  the  depth  was  within  ten  ireet  of  the  ordered  depth.  The  model 
was  Initialized  on  course  000  degrees  true  at  an  Initial  speed  of 
twenty  knots. 

The  maneuver  was  successfully  completed.  The  model  remained 
steady  at  320  -1  degrees  true  and  the  final  depth  was  within  one 
foot  of  the  ordered  depth.  Both  the  pitch  and  roll  angles  were  near 
one  degree  and  were  decreasing  in  magnitude.  All  of  the  control 
surfaces  were  at  zero  angle  of  deflection.  It  should  be  noted  that 
the  new  course  had  been  achieved  during  the  first  sixty  seconds  and 
that  the  course  was  maintained  by  the  model  while  the  proper  depth 
and  attitude  were  being  obtained. 

IV.  USE  AS  A  DESIGN  TOOL 

There  are  nunerous  design  tasks  that  could  be  used  to  demonstrate 
the  simulation  model;  a  single  example  Is  sufficient  for  this  demon* 
stratlon.  Suppose  a  designer  wanted  to  know  the  dynamic  effects  of 
Increasing  the  rudder  rate  In  a  cum.  The  designer  would  elect  to 
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use  a  simulation  model.  Since  he  would  not  want  the  dive  planes  to 
interfere  with  the  analysis,  he  would  set  NOOIFF,  DCRIT,  and  N0P1CH 
a4,  large  values  so  that  the  dive  planes  would  not  move.  Then  he 
would  set  RRATE,  the  rudder  rate,  at  the  value  he  d&slred  to  test 
and  order  the  model  to  come  to  a  new  course.  For  this  run,  let's 
as suae  that  the  designer  performed  the  left  hand  turn  from  000 
degrees  true  to  320  degrees  true.  He  used  a  rudder  rate  of  two 
degrees/second.  He  then  caused  the  model  to  perform  the  same 
manuever  again  with  a  new  rudder  rate  of  four  degrees/second.  With 
the  output  from  these  two  maneuvers  he  could  easily  find  the  time 
required  for  the  turn,  the  advance  and  transfer  of  the  submarine, 
and  the  roll  and  pitch  angles  as  a  function  of  time.  The  designer 
could  use  the  information  for  whatever  analysis  he  had  In  mind.  The 
model  could  be  run  again  at  new  rudder  rates  or  the  dive  planes  could 
be  brought  into  play  or  virtually  any  other  maneuver  could  be 
simulated.  This  model  can  also  do  snap  roll  analysis  [4]. 

The  cost  of  running  the  simulation  model  Is  so  low  that  It  can 
be  used  on  a  dally  basis  If  desired.  For  the  test  turns  mentioned 
above,  the  cost  was  less  than  five  dollars,  which  included  reading 
the  cards,  compilation,  execution,  and  900  lines  of  output.  A 
designer  who  used  the  model  regularly  could  have  the  model  as  an  on 
Una  dataset  which  would  greatly  reduce  the  cost  of  using  the  model. 

IV. 6  CONCLUSIONS  ANO  RSCOWEN  OAT  IONS 

This  simulation  model  does  accurately  simulate  the  six  degree 
of  freedom  motion  of  a  submarine  for  moderate  maneuvers.  The  cost  of 
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operating  the  model  Is  very  Iom,  which  will  «11m  frequent  usage. 

The  model  Is  designed  to  permit  a  great  deal  of  flexibility  In  the 
application  of  the  model.  For  those  designers  who  use  the  simulation 
model.  It  should  be  a  useful  design  tool. 

If  work  were  continued  on  this  simulation  model.  It  Is  recam- 
mended  that  some  time  be  spent  In  Improving  the  appendage  control 
subroutines.  Further  application  of  control  theory  would  be  helpful 
with  the  appendage  subroutines.  The  data  Input  could  be  organized 
more  efficiently  to  remove  some  of  the  opportunity  for  error.  It 
could  be  very  useful  to  have  a  plotting  routine  In  the  model  to 
visually  display  the  Information. 
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A.1  NOTATION 


Symbol  Dtmeneionlees  form 


Definition 


B 


Buoyancy  (ore*.  ponUvg  upward 


CB 

CG 


t 


I  •  .  »* 

“  tot 

VKP- 
s  •  -4=V 


i  ■  .A. 

r» 


Cantor  of  buoyancy  of  submar-ne 
Cantar  of  min  of  submarine 

Moment  of  inertia  of  submarine  about  x  axis 

Mamant  of  mama  of  submarine  about  y  axia 

Uomaat  of  tnartia  of  submarine  about  a  axta 

Product  of  martta  about  xy  axta 
Product  of  martta  about  yt  axaa 


t 


ax 


K 


% 

K 

P 


r. 


Pd 


ax 


*ti 

txT 


K'  • 


K  '  e 


i'bt5  u* 

K. 


TTTI 


"  "  i!»l 
^  tel  w1 

K  •  •  ■*? 


P  *pt*u 


K.-  .-5L 

>  t«*f 


K 

Pd  iti‘ 


Product  of  tnartta  about  ax  axaa 


Hydrodynamic  moinont  componant  about  x 
aai  trolling  moment) 

Rolling  moment  »h«i  body  angle  in.  3)  and 
control  surface  anklet  are  tero 

Coefficient  used  m  representing  K«  aa  a 
function  of  (rj-  l) 

Tirat  order  coefficient  uaed  in  representing 
Kan  (unction  of  p 

Coefficient  used  in  repreeenting  K  as  a  function 

of  p 


Second  order  woeificient  used  in  representing 
K  at  a  function  ot  p 

Coefficient  jv:d  ir  representing  K  aa  a  function 
of  the  product  pq 


V 

K*r  *  i*1" 

Coefficient  ue»d  in  representing  K  as  a 
(unction  of  the  piMuci 

*r 

First  order  coefficient  used  tn  representing 
Km  function  of  r 

*'m& 

Coefficient  used  in  representing  K  os  * 

(unction  of  f 

*v 

*  .  K* 

*  ’U*u 

First  order  coefficient  used  in  representing 

K  ss  a  function  of  * 

S 

„  ,  K* 

1C*  ■  7 - 

’  io** 

Coefficient  used  in  representing  X  as  a 
function  of  » 

K*l*l 

„  Kvl,l 

•  *1*1  *  F7Z’- 

Second  order  coefficient  used  in  representing 

K  as  a  funcron  of  » 

**, 

K  • .  Jba. 

**  iet* 

Coefficient  used  in  representing  K  as  •  function 
of  ihs  product  rq 

Kra. 

*  .  *n~ 

K  »  r— t‘ 

"•  ipT* 

Coefficient  used  in  representing  Kit  a  function 
of  the  product  « 

K-p 

K  • .  J5a 

*»  »»4‘ 

Coefficient  used  in  representing  K  is  a  function 
of  ihr  product  »p 

*wr 

■W-Bff- 

Coefficient  used  in  representing  Kiel  function 
of  the  product  »r 

*»r 

k  •  . .  *$r.  | 

»r  *p4JU* 

First  order  coefficient  used  in  representing 
Kiel  function  sf  if 

4 

4*  •  l 

Overall  length  of  submarine 

m 

«.•  .  — 2U 

Mass  of  submarine,  including  meter  in  free* 
flooding  spacee 

M 

“'•CTF 

Hydrodynamic  moment  component  about  y  aula 
(pitching  moment) 

*• 

Um 

*•*  ’^Tu* 

Pitching  moment  »mn  Sod v  angles  (<».  S)  end 
control  surface  ingles  ire  tero 

*PP 

M  -.r^ 

PP  *®4‘ 

Second  ordor  coefficient  usd  .n  representing 

M  as  a  function  of  p.  First  orosr  coefficient  is 

sera. 

M 

* 

M  ■  - 
*  |®4*U 

First  order  coefficient  used  in  representing 

Id  as  a  i unction  of  q 

W_ 

Vt 

m  •  .  JSia- 

First  order  coetflcient  used  m  representing 

M  as  s  (unction  ol  o -  ) 
u 

M. 

d 

Mi 

Ml*'  •  v- ~7 
*  t®4» 

Coefficient  used  hi  repr  tenting  Mass 
function  of  q 
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Second  order  coeKicionl  u*ed  \n  rep r esc mi n| 

M  4*  4  (unction  o<  q 

Coefficient  used  in  representing  »•  4 

function  q 

CMdicii't  used  ir.  representing  M  as  * 
function  of  tho  product  rp 

Second  nrd' •*  cocf'icient  used  in  representing 
M  4i  a  function  of  r  first  order  cooificicat 

it  tors 

Coeificienl  used  in  reprcienttng  Mae* 
function  oi  the  product  vp 

Coefficient  need  m  representing  M  at  a 
function  ot  .he  product  vr 

Second  order  coefficient  need  in  representing 
Mata  function  of  * 

rtrat  order  cocificient  used  in  repretenttng 
M  at  a  function  of  » 

rtrtt  order  coefficient  used  in  repretenttng 
Mw  «•  t  function  of  (n-  l) 

Coefficient  used  in  representing  M  as  a  function 

oi  » 

First  order  coefficient  used  m  representing  M 
at  a  function  of  equal  to  icro  for  lymmetncal 
function 

Coefficient  used  in  representing  M_  at  e  function 

oi  • 

Second  order  coefficient  used  n  representing 
Matt  function  of  » 


rtrtt  order  coefficien:  reed  m  representing 
M  ,  ,  a>  a  function  oi  in-i) 

Second  order  coefficient  used  in  representing 
M  as  a  function  of  » .  equal  to  tero  for  iftn- 
metrical  function 

First  order  coefficient  ited  in  representing 
M  at  a  (unction  of 

First  order  coefficient  used  in  representing 
N4  ss  s  function  of  5  , 

First  order  Cieffictert  used  in  representing 
•Vis,  u  i-mcion  ot  id-  l) 
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N 

^  ** 

"  *  rjrqp 

M>  drodvn.erme  moment  comoorifni  about  a 
arts  (yawing  moment) 

N. 

N*  'T&'u' 

Yautng  moment  uoen  Sort y  injlee  [ft.  3)  and 
control  turftce  -inglet  are  tero 

.X 

? 

V  ’tS^r 

rirat  order  coefficient  utod  in  rtpreaenling  N 
aa  a  (unction  o(  p 

N.«.  A- 

l»  to** 

Coefficient  need  in  repre  tenting  N  at  a  (unction 
of  £ 

N 

PX 

N  *  •  -^23- 
n  ***‘ 

Coefficient  need  in  repretenting  X  at  a  (unction 
of  the  product  pq 

N 

jm  A 

K  '  *  » - 

^  to** 

Coefficient  u»cd  m  repreteniing  N  at  a  function 

of  the  product  qr 

M 

r 

N 

Nr'f^U 

Tirtt  order  .ouKicient  uted  in  repreaenttng  X 
aa  a  (unction  o(  r 

Nn» 

m  • .  J^n. 

rn  tO**y 

Tlrtt  order  .  ^efficient  uted  in  repretenting 

Xf  at  a  funci.on  of  f *y-  1 1 

m$r 

Coefficient  uted  •.  i  repretenting  Mata  (unction 
of  r 

Vr|rt 

\...  -55^ 

Second  *•--*.-?  coefficient  uted  in  repretenting 

N  at  a  function  of  r 

*!rl»r 

Nl ri «r 

X.  .  •  e  - L_ll 

rl4r  jofU 

Coeffleie n*  .»ed  .i.  repretenting  it  a 

fuiwuot  of  r 

X 

* 

N 

X  •  .  , - S_ 

"  to*  U 

Tlrtt  nrdtr  cceffic;  tr.t  uted  in  repretenting  .V 
at  a  function  of  v 

M 

x__ 

x  •  .  — -23L 
•*»  to*  u 

n-tt  order  coefficient  uted  in  repretenting  Nv 
at  a  function  of  ,q-l) 

H. 

» 

’  to** 

Coefficient  u»t  d  .n  repretenting  X  tt  a 
function  of  * 

N 

x  -.-Sa. 

^  |o** 

Coefficient  ut*d  .n  repretenting  V  tt  a  function 
of  the  product  vq 

K'.l, 

M,  ■  . 

!*lr  ,*• 

Coefficient  uted  in  repretenting  Vr  at  a 
function  of  » 

M»l*l 

x  •  - 

''<*1  to*3^ 

Second  order  c  ocfficunt  uted  m  repretenting 

N  at  a  function  of  v 

N»iv(n 

n  ,  ,  • . 

"I*1*?  |S4' 

FlrJt  order  coetficient  uted  in  repretenting 

X  vl  »  t  t  nine*. ion  of  .-J-  1 ) 
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Coefficient  u>ed  in  ftjrnentint  N  it  •  (unction 
of  the  product  « 


„  . 

N  •  f  "■  'gl 


Coefficient  used  >n  representing  N  He  a  function 
of  the  product  «p 


Wf 

f  • 

.  ,  "&r  , 

dr 

f 

N)n, 

bn) 

7 


Coefficient  uacd  in  representing  N  ae  a  function 
of  the  product  *  r 

Tiret  order  coefficient  used  in  representing  N 
ae  a  function  of  ir 

first  order  c  f.cient  used  in  representing 
N jr  aa  a  function  of  !r)-  l) 


Angular  velocity  component  about  \  i*ie 
relative  to  flu  i  (toll) 


Angular  ae  ce  .ration  component  about  v  rata 
relative  to  fluid 


Angular  velocity  component  about  y  axle  relative 
to  fluid  (pitch) 


A  ngule  r  eccnioretitn  ttmponin;  about  y  me 
relative  to  fluid 


Angular  ve.x.tv  component  about  t  sxta 
relative  to  fluid  :  yaw) 


Angular  act  elvratinn  component  about  i  i*i» 
relative  to  fluid 


Linear  velocity  of  crigm  of  body  aaee  relative 
to  fluid 


Component  of  V  in  diftction  M  th#  i  ixli 


Tim#  nti  of  of  o  in  'ftroetton  of  th# 

a  axil 


Command  tp*.**4,  •  i  c  *  ^  v  m |  j*  of  imH  ip««d 
compontnt  n  [jt  a  t'llr r  rpm  ih«A 

'i.  !'  mr*  .ontroi  turUcn  -ingle# 
*uro.  S.|n  chAnfot  %uh  propeller  ’’evertAl 

Component  of  V  .  n  4tr#cuon  of  the  y  axis 


Tim#  r»i#  of  chetnt**  of  *  to  'Mnrtion  of  th# 

r  «a;# 
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Component  <?(  U  in  dirictton  of  the  i  axis 


Tim*  rate  of  change  o(  »  tn  direction  of  the 
i  udi 


W 


W 

toPtF 


Weight,  including  water  in  free  flooding  apace* 


Longitudinal  body  ini,  alto  the  coordinate  c(  a 
point  relative  to  the  origin  of  body  axaa 


Tb#  a  coordinate  of  CB 


The  a  coordinate  of  CC 


A  coordinate  of  the  displacement  of  CC  relative 
to  the  origin  of  a  set  of  fixed  axes 


X' 


Hydrodynamic  force  component  along  x  axis 
[longitudinal,  or  axial,  , oriel 


Second  order  coefficient  used  in  representing 
X  a*  a  (unction  of  g.  First  order  coefficient 
is  sere 

Coefficient  used  in  representing  X  ae  a  function 
of  the  product  rp 


X.' 

u 


Second  order  coefficient  used  in  representing 
X  aa  a  function  of  r  first  order  coefficient  te 

aero 

Coefficient  used  in  representing  X  ae  a  function 
of  u 


X 


uu 


*sm 


Second  order  coefficicni  used  in  representing 
X  at  a  function  of  u  in  the  non- propelled  case, 
first  order  coefficient  is  sere 


Coefficient  used  in  re;* -e tenting  X  ae  a  function 
of  the  product  »r 


Second  order  coeffient  used  io  representing  X 
aa  a  function  of  ».  First  order  coctlicte  n  it  zero 


First  ord»r  coefficient  used  in  representing  Xvv 
ae  a  (unction  of  (rj-  1) 


Coefficient  used  in  representing  X  at  *  fu  ictton 
of  th*  preduct  •<! 
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X.w 

X  •  e  ,^L-L 

erer 

x  •  -fel 
wwn  tuf.1 

X6bAb 

w  .  .  xAb4b 

4b4b  gpA,U1 

X4rdr 

X.  *  - 

4r4r  fpt^1 

X$r4nj 

v  ,  _  X4r4rrt 

Ar4rn 

XSs4e 

• 

y  .  XAsAs 

4s4a  | p4 ; U 1 

X.  * 

4  s  A  sty 

y  .  _  *4e4e»j 

4s4stj  ipg,*!.'1 

r 

y*  a  X 

t 

r3 

V  ■  T*- 

fC 

rG 

yG  •  T 

ro 

r0 

y  ■ 

Y 

r 

'  frpC'U* 

Y.  e  * 

r 

p 

Yp 

Tp  *  ice*  u 

Yp 

Y.-  •  r^i~ 

P  |P4* 

V,' 

V  .  _  Tp 1 »  t 

pip1  gp4« 

Second  order  coefficient  used  in  representing 
X  a*  *  function  of  *  .  '  First  order  coefficient  is 

tero 

First  order  coefficient  used  in  representing 
as  a  function  ct  (n -  l) 


Second  ord"  coefficient  used  in  representing  X 
as  a  function  of  First  ordor  coefficient 

Is  tero 


Second  order  coefficient  used  in  representing 
X  as  A  function  of  4r.  First  order  coefficient  is 

s*s>rs  * 


First  order  :o*ffi>ient  used  in  representing 
Xjr4r  **  4  ncuori  (n*  1) 


Second  order  coefticient  used  in  representing  X 
as  a  function  of  i,.  First  order  coefficient  is 

sero 

First  order  coefficient  -used  .n  representing 
*5sjs  **  4  function  of 

Lateral  bod v  ixii;  Also  the  coordinAte  of  a 
point  retail-*  ia  the  origin  of  bed y  islt 

The  y  coordinate  <-f  CO 


The  y  coordinate  of  CC 

A  coordinate  of  th  i  displacement  of  CC  relative 
to  the  origin  of  s  set  oi  fixed  sxes 

Hydrodynamic  forte  component  Along  y  axis 
(Lateral  force) 

Lateral  force  whan  body  angles  (a,  3)  and  control 
surface  angles  ire  irro 

First  ordor  coefficient  is-td  in  representing 
Y  ss  s  func'ton  o;  p 

Coefficient  used  in  representing  Y  as  a  function 

of  p 

Second  order  coefficient  used  in  rep  re  sentir.g 
f  til  function  of  p 
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V 

Y5»5 

.  Jju. 

*ct« 

Coefficient  ueed  tn  rcpreeentinj  Y  \»  a  (unction 
o(  the  product  p<\ 

r*r 

.  J!ii 

«pf 

Coefficient  uted  in  repreeentinj  Y  ae  a  function 
of  tUa  product  qr 

Yr 

Yr'  * 

rsiH 

Flr»l  ordor  coefficient  uted  in  repretoniin|  Y 

4*  4  function  of  r 

rm 

*f  * 

*  »~n 

‘FTFu 

Flrtt  ordor  coefficient  ueed  tn  representing 

Yr  4»  4  function  of  (rj-  l) 

Tf 

V  ’ 

** 

tpt‘ 

Coefficient  ueed  in  repreaenttng  Y  ee  4  function 
oi  i 

Ylrlir 

-  ,  T’H»r 

Ylr!J»  '  ilV'T 

CooffUiont  u»*d  in  reproeenting  Y  4t  4 
function  of  r  r 

Y 

Y  '  * 

jr. 

Flrtt  ordor  coefficient  ueed  in  repreeenting 

Y  44  4  function  oi  v 

v 

* 

SOT 

% 

V 

.  J^L 
ta4'« 

Flrtt  ordor  coefficient  uttd  .n  repreaenung 

Yv  4i  4  function  of  i.,-  i) 

n 

T .  *  • 

* 

*v 

Coofftciont  utod  in  eepreeentlng  Y  u  4 
function  of  v 

% 

T  ’ 
*4 

Coefficient  ueed  tn  rtprtitnnn|  Y  to  4  function 
of  tHe  product  «( 

T*'r( 

w  ”» 1  r 1 

T  ,  e  -r  ■■  -* 

*i  H  |o<. 

Cooffloont  utod  tn  r»prti«nltn|  .  v  aa  a  function 
of  r 

Y»lvf 

Ye<e 

.  T.lrl 

•  "  ,0<r 

Second  order  coefficient  ueed  in  repreeenttng 

Y  it  i  function  of  v 

Yvl V 1 ») 

Y,i, 

^  vUltf 

'•»  'TIT1-’ 

Flrtt  order  coefficient  ueed  tr  rtpreeenling 

Y  y  i  v  |  4t  e  function  of  rj-  1) 

r 

Wf 

Y~ 

*  *Pt‘ 

Coefficient  ueed  ,n  rroreeennag  Y  it  i 
tuncton  of  the  product  v» 

Y 

T»F 

T_ 

Coafflciant  ueed  in  reproeenting  Y  ee  a 
(unction  of  the  product  »p 

*r 

T 

m  . 

«  I*LL 

tPt 

Coeifictent  ueed  tn  repretennng  Y  ee  4 
function  oi  tiie  product  » r 

T»r 

Ylr 

•  J^r 

|Ot*U‘ 

riret  orcer  coefficient  ueed  in  reproeenting 

Y  at  a  (unction  of  Jr 

T»r» I 

Y»rr, 

•  .  -.Jja. 

»pt*u* 

Firet  order  c-ocffictcr.:  ueed  in  repruaetiting 

Yjf  ee  e  iunciion  oi  In-  U 
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s 

L 

£ 


zwi»i 


*w|w|tj 


_  ^wlwl 

Z-M  "  tZF~~ 

m 

z»'»!n  ”1^ 


'M 

z  • 

are 

■nrr_ 

*  • 

*»k 

ftb 

4b 

*  TpPu7 

r.  • 

•  Zl* 

ft# 

fte 

I.  1 

*fteie 

**n 

*«n 

urn? 

* 


T 


■ 


S*cORd  order  coefficient  used  in  representing 
Z  41  l  function  of  e 

rirst  order  coefficient  used  in  representing 

Z  .  ,  as  a  function  of  (n-l) 

*|w! 

Second  order  coefficient  used  in  representing  • 
Z  is  i  function  of  w;  equal  to  tero  for  sym- 
metricsl  function 

rirst  order  coefficient  used  in  representing  Z 
ss  s  function  of  4h 

First  order  coefficient  used  in  representing 
Z  aa  a  function  of 

First  order  coefficient  used  In  representing 
Zj#  ss  k  function  of  (q-l) 

Angle  of  stuck 

Aegis  of  drift 

Deflection  of  bowplane  or  teilplsne 
Deflection  of  rudder 
Deflection  of  sternplsne 

Ur 

The  ratio  -~- 
Angle  of  pitch 
Angle  of  yaw 
Angle  of  roll 

Sets  of  constants  used  in  the  reprs sentation  of 
propeller  ibrust  in  tne  ssisl  equation 


A. 2  EQUATIONS  OF  MOTION 
AXIAL  FORCE 


m^u  •  vr  ♦  wq  -  xc  (qa  ♦  ra)  +  Yq  (PR  •  *)  +  %q  Ipr  *  q)]  » 

♦  —  f  Tx  '  qa  +  X  ’  ra  +  X  '  rp  1 

1  L  qq  ^  rr  rp  r  J 

♦  ~  i*  ["*■  '  u  ♦  X  *  vr  +  X  '  wq  1 

2  L  u  vr  wq  ^  J 

♦  —  la  [x  '  aa  +  X  *'  v2  ♦  X  ’  w*  1 

2  l-  uu  vv  vvw  J 

♦  —  i*ua  fx.  *  '  5ra  +  X.  •  6sa  ♦  x..  '•  5ba 

2  L  6r  6r  6s  6s  5b  6b 

+  1  pta  «a  +  b.  auc  +  c.  uca  j 


-  (W  -  B)  sin  9 


7-  ta  T  X  va  *  X  ’  wa  ♦  xt  ,  '  6  a  aa 

2  L  vvTj  #wTj  Sr  Sr  q  r 


NORMAL  FORCE 


,[w  .  u,  ♦  vp  -  *c  <P*  ♦  1a>  *  *c  1'P  -  <>  *  *G  tr’  *  **»  j  ■ 

*  S 

♦A  '*  K' <  *  V 

♦  ±  s  [zw'  «  ♦  zvr'  «  *  V  *P  ] 

♦f  «*  K'  «!  *  z„l4.'  “'^5*  *  w  ^  '"J 

*j  t*  [z.'  a*  ♦  zw'  aw  *  ZW|W|'  "'(V* 

♦■f  *’  Cz|w!'  alw|  *  z-~'  |w  ty’  + 

+f  ia  [z^1  va  ♦  Z^1  u*  6s  ♦  Z6b'  ua  6b]^ 

♦  (W  -  B)  cos  0  cos  * 

+  f V 

*4  t*r*.  '  «»  +  Z„l„ln'  “'ll*’  +  ZWV’j(,M1 


ROLLING  MOMENT 


♦  m 


(I.  -  I J  q r  -  (r  ♦  pq)  Ixa  ♦  (r3  -  q3)  ♦  (pr  -  q)  I 


*  y 


y* 


xy 


|^yG  (w  -  uq  ♦  vp  )  -  (v  -  wp 


ar)  ] 


-L  XBI  K  '  p  ♦  K-  '  r  ♦  K  '  qr  ♦  K  '  pq  +  K  ,'p|p 
2  Lpr  r  qr  pqr^  pipI  r,r 


pq 


pipI 


-  X^K  '  up  ♦  K  '  ar  ♦  K  '  v  ! 
2  1  P  r  v 


2  ~  L  p 

♦—  l4  .rK  '  vq  *  K  '  wp  ♦  K  •  wr  "J 
j  L  vq  ^  wp  r  wr  J 


♦-2-  X3  [v  u3  ♦  Kv'  uv  ♦  Kv|vj'  V  |(v3  ♦  w3,*|l 


♦J.  x3  I  K  'vw  ♦  K  ’  u3  Sri 
2  i.  vw  Jr  j 


♦  ( yr  W  -  yB  3)  co*  6  co*  ^  -  (*GW  -  *bB)  co*  8  *in  ♦ 


B 


+  f4*  K*T7’  u*  M-U 
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.3  AXIS  TRANSFORMATIONS  [2] 

1)  A  Crans formation  of  an  axes  system  takes  a  quantity  descried 
In  one  frame  of  reference  and  transforms  It  Into  another  frame  of 
reference  such  that  If  we  measured  the  same  quantity  In  the  second 
fram  of  reference  the  transformed  quantity  and  the  measured  quantity 
would  be  Identical . 

2)  Transforms  between  frames  are  needed  In  the  study  of  the  motions 
of  ocean  venlcles  because  the  equations  of  motion  for  such  a  vehicle 
are  most  easily  derived  In  the  Inertial  frame  attached  to  the  earth 

(x  ,  yQ  ,  zQ)  frm m,  while  the  forces  acting  on  the  vehicle  are  most 
easily  evaluated  In  the  frame  attached  to  the  vehicle  (x,  y,  z). 
rtence,  we  ultimately  desire  to  transform  the  equations  of  motion 
frm  the  Inertial  frame  Into  the  non-inertlal  frame  fixed  In  the 
vehicle. 

3)  If  V0  Is  some  vector  measure  in  the  xQ  ,  yQ  ,  zQ  frame  and  V 
some  vector  measured  In  the  x,  y,  z  frame  wnlch  Is  only  cnanged 

In  orientation  then: 

V  •  T  (  0 ,  * ,  *  )  wnere  T  (*,-*.  ®  )  -  the  transform. 

Where 

j"cas*  cose 

!  -sin*  cos®  ♦  sin®  sin*  cose 

i 

!  sin®  slnv  *  cos®  cos*  sin* 

-sin  *  | 

sin®  cos-*  ! 

COS®  COS'* 


T  (e,  -*,  *  ) 


cos*  sine 

cos®  cose  ♦  sin®  sin*  sine 
-sin®  cose  ♦  cos®  sin*  sine 
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8.1  LIST 
VARIABLE 

A 

AA 

A1  ,A2 

AI.8I.CI 

8 

BORATE 

BOMNAX 

COURSE 

OECAIT 

oai 


APPENDIX  8 


OF  VARIABLES 


MEANING 


Six  by  six  matrix  containing  coefficients 
Of  UQT,  VDT,  WDT,  POT,  QDT,  ROT. 

Six  by  six  matrix  set  equal  to  matrix  A 
before  each  call  to  FUNC.  The  values  of 
AA  are  lost  In  the  matrix  reduction  performed 
by  LEQT1F. 

Limits  used  for  selecting  the  proper 
propeller  thrust. 

Set  of  constants  representing  the  propeller 
thrust  In  tne  X-equation. 

Ship's  Buoyancy. 

Average  bowplane  rate. 

Maxlwjn  ordered  oowplane  deflection. 

New  course  for  the  snip. 

Value  of  depth  error  when  dive  planes  are 
returned  to  zero  deflection. 

Calculated  deflection  of  the  bowplane. 


DEL  3 


Calculated  defection  of  the  stem  plane. 
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0EL4 

OELB 

oaeo 

D£LOi 

DELR 

QELRO 

DELS 

OELSD 

oar 

OIFF 


IC.NT 


Calculated  deflection  of  the  rudder. 

Actual  bowpUne  deflection  at  time  t. 

Actual  bowplane  deflection  (units  changed 
for  output). 

Amount  the  original  GM  Is  to  be  changed. 

Actual  rudder  deflection  at  time  t. 

Actual  rudder  deflection  (units  changed 
for  output). 

Actual  stem  plane  deflection  at  time  t. 

Actual  stem  plane  deflection  (units 
changed  for  output). 

Time  Increment  used  In  Iteration. 

Difference  oet*een  present  depth  and 
ordered  depth. 

Six  by  one  matrix  containing  solution  to 
right  nand  side  of  equations  of  motion. 

Counter  used  to  coieit  the  nietoer  of 
1 1*  rations 


ID  Forty-character  alphanumeric  heading. 

INDEX  tf  greeter  than  zero,  read  new  submarine 

coefficients.  If  less  than  or  equal  to 
zero,  run  same  submarine  for  new  initial 


conditions. 
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IX,1Y,IZ,IXY , IXZ.IYZ 
KCOEFF 

L 

LEQT1F 

H 

MAXANG 

HCOEFF 

NCOEFF 

MOOIFF 

MOPtCH 

OOEPTH 

OMCRS 

P,Q.» 

POT.QOT.aOT 

POTO.QOTO.RDTQ 


Moments  of  Inertia. 

K  -  equation  coefficients: 

Ship's  overall  length. 

Matrix  reduction  subroutine  from  IHSlIB  [3] 
Ship's  mass. 

Maxlaua  ordered  dive/ ascent  angle. 

M  -  equation  coefficients. 

N  -  equation  coefficients. 

Acceptable  error  range  around  ordered  depth 

Acceptable  error  range  around  zero  pitch 
angle. 

Ordered  depth. 

Acceptable  error  range  around  ordered 
course. 

Angular  velocity  about  the  x,  y,  and  z 
axes  respectively . 

Angular  acceleration  about  Lfte  x.  y,  and  2 
axes  respect! vely. 

Angular  acceleration  about  tne  x,  y,  and  z 


PHI 

PHIO 


PO.QO.RO 


PS  I 

PSIO 


RHO 

RRATt 

RUCAHT 

STERAT 


STtRHX 


STOW 

T 

Ti  .  .  .  na 
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Axes,  respectively  (with  units  changed  for 
output). 

Angle  of  roll . 

Angle  of  roll  (with  units  changed  for 
output). 

Angular  velocity  about  the  x,  y,  and  z 
axes  respectively  (with  units  changed  for 
output ). 

Angle  of  yaw. 

Angle  o'  yaw  (with  units  changed  for 
output) . 

Sea  water  density. 

Average  rudder  rate. 

Maxima  ordered  rudder  deflection. 

Average  stem  plane  rate. 

Maximal  ordered  stem  plane  deflection. 
Storage  array  for  naif  of  output  data. 
Present  time. 

Time  lag  signals. 


THETA 


Angle  of  pitch. 
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THETAO 

TU68 

TLAGR 

TLA6S 

u.v.w 

UOT.VDT.WOT 

uo 

uoo,vo.*g 

*T 

X.Y.Z 

Xfl,Y3,iB 

XCOEFc 

XDT,YUT,Z£T 

XG.YG.ZG 


Angle  of  pitch  (with  units  changed  for 
output). 

Tine  lag  for  Powplane  control  system. 

Time  lag  for  the  rudder  control  system. 

Time  lag  for  stem  plane  control  system. 

Forward,  lateral,  and  vertical  velocities 
nespectl vely. 

Forward,  lateral,  and  vertical  accelerations 
respectively  (wltn  units  changed  for  output) 

initial  forward  velocity. 

Forward,  lateral,  and  vertical  velocities 
respectively  (with  units  changed  for  output) 

Ship's  weight. 

Coordinate  laoels  of  the  fixed  (xQ,  yQ, 
zQ)  coordinate  system. 

The  x,  y,  z  position  of  the  canter  of 
ouoyancy. 

*  -  equation  coefficients 

Velocities  in  the  fixed  coordinate  syltem 
along  the  xQ,  y  ,  z^  sixes  respectively. 

Hue  x,  y,  z  position  of  the  center  of 


gravi ty. 


YCOEFF 

Y  - 

ZCOEFF 

Z  - 

equation  coefficients 
equation  coefficients 
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8.2  COMPUTER  PROGRAM  LISTING 


kLAU(b ,10  (ZCOCFFU  ,1-1,22 

RLAI)(b,10  ICCOtFF  1  .1-1,17 

kLAO(b  ,10  MCOEFF(I).I-l .23 

klAD(b,lO  HCOEFF( I  ,1-1,22 


KlAOfb.lO 

k£Au(t>,10 


klAU  UtPTH  CHAMGt  INFORMATION 


IS)*  ,6X,‘  (KTS)'  ,4X,'  (XTS/itC)*  ,2X,'  (KTS/SIC)  *  ,2X ,  ‘  (KTS/SEC) 1 ,3X .  (0 
2t6)'  ,7X,'(Ut6)' ) 


AC 2  ,4}»-(M*Z6*YC0€ff  (2  )*L) 
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T-O.OUl 


64 


o  o  o 


THLTA- ML  IA/67 . 2ybd 
PSl-PSI/b/^btJ 
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i 


k-k*utLT*Ror 


x  — •  — 

2  x  x 

t—  X  X 

>-  i/i  yt 

✓»  ©  o 

Ui  XX 

►-  •  ♦ 

5  ~ 

M  l/»  V/1 

g  2^ 


□  i/ti/i 

3  o 

a  •  « 


« i/i  &. 
i/i  o  — ■ 
3  «  /> 

♦  —  O 

-5  3 


♦  a.  ♦  x 


ox  —  «  ^ 

X^  3?  ~ 


l/l 

C  /I  — < 

—a 

vl 

^ 

on  >— 

u. 

«kJ  '  _ i  1*4 

O.  JZ 

3 

X  1/1  3  X 
•-Oil- 

X  3  — *•—• 

3 

</i  35  - — — j 

o*  -oa''i/i^ 

ST**  —  5?? 

■— *  ►-  >—  — »  <  i/l ' — -  •*“** 
•  oJ  uJ  •  X  * 

§g£3gg3g 

*5  5i  3 

/i  Vi  — *  >— 

3  0  i£— -  — 

a  a  •  a  —  <  a  ^ 

3  3  ■ — '333 

»«X/1^-«X'>40 
—  M  *.»■  na  — 


—  i/>  <. 

i/t  a. 

a.  ♦  x 

J  l/l  O 
/IX  3 

5 

sr  5 

/I  i/l  3 


a  ■  3 

•Ni  .■Nl  3  __ 

i/l  i/l  r»  3 

XX  X 

o  u  />  «C 


2*  2*2*7 

*>  *  -k. 

^ 

a  a  a  — 


—  x  •  • 

X  i  3  3 

X  -J  ■  • 

Cl  1^  3 

i  l/l  /I  X  _J  o 

I  sa  —  3  •  • 

:<^3  —  — 
i  3  •—  O  i/l  vl 

a  •  a  X  X 

>  —  I~1  —  3  3 

i  x  a  /I  — • — • 
,  x  —  x  3 
:  <  x  o  —  — 


5  Zi 

X  ~ 

X  X* 

— J 

3  I 


*.£5* 

X  uJ  3  >— 
3  3  ►-  < 

3  *r-  — 

ari-,3 


■/■'  3  ►—  — 
X.  —  X 

o  a:i  x  o 

oii  • 


J  3  O 


3  3  3  < 


uuu 


0  0  3  3 


COHVEKT  AKGLtS  1M  RADIANS  TO  OiCM.ES  AMO  FT/SEC  TO  Wfc)TS  FOR  OUTPUT 


UHl  It.  (6  ,bO)  T,  X  .1  ,U00 .  VO  ,tt).UOTO.  VUTO  .WOTO.DtLSO  ,D£LBO 
it)  FWJHAI(3(U,010.3).till.4e8O*0l0-l)./) 


SUbKOUHNt  FUNC 


CJI-CK2) 
GO  10  40 


imyfJMW  IP  l|«JUiJ|p!|^JLl 


3*  */l  & 

sis 

3?5 

-*4  *—  *— 

♦  f\l  — 
y^>  ^ 

-I  U.  b  3* 

r>  o  —  3 
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u.  ac  a»  — -a — * 

,U  «  ~— '  —  bN  fvi 

8? x<42 

M  — ■«<  ►—  a  O' 
— »  _  o  Mi  •  ■ 

•  . —  a  wON 


.—  Ml  ^  «  — — — — 

—  O  —  a  ac  — ■ — 

~—  U  ■ — ■  ±3  —  i  — • 

♦  x  u.  —  x  —  a. 

—  ♦  —  a  a.  a.  • 

—  a“  <5>  —  £  O' 

33  S  X  —  O  -N  • 

3  —  *  —  a  »<  o 

a  O  O  •*•  •  —  •— 

•  —  •  mI  — « 

a.  — ■  >  3  «  ♦  o 

a  u.  a  i_j  <  ■ — -  'n 

—  —  — ■  Sna(  I 

mI  i">  ♦  Mi  •  ^ 

—  3  —  ac  x  crry 

m.  W.  —  Mi  — -—  3 
a)  ♦  u.3^>  i 

o  o  w  «  o  —  a. 
o  •  o  — O  •  • 

«;»uNaN> 

♦  •  2  a  «  —  — 
O’  —  ♦  •  — — '  « 

•  er*  —  o  a  •— -  * J 

a.  • —  cm  ■ —  «  i  »■ 

•  u.  «  «  a  - — ■■ — 

- — ■  —  •  — >»  —  a 
«  w  ^  -O  1  — 

—  o - t  at  — 

u.  O  a  —  a  —  cm 
u.  2  —Mb  •  X  a 
lm  *  (m  m  '  J  a.  • 
OS-mk-J 
L_1  •  — 'O 

2  o  m.  o  — —  — 

*  a  m.  2  a  v-">  - — • 
0£  -*bj  »  —  O  • 

a  —  O  3  — «  X 
fyw  <_)  •  ^  —  -w' 
t  u.2><  at  — 

—  —  —  a  a  C  i 
O  Mi  a  —  j  h-  CC 
— -O' — ■  u/1  «  mJ  a 

bUN-OSl 

—  2  a  -*Ha 
^  *  a  —  2  —  x 

o  a.  _i  —  —  v;  x 

SJ  a  "v  £3  — 
2  3  -  3  n  y  < 

«  a  —  vj  — -.O  — 
-m— .  ^  —  a  *N 

a  a  ♦  *  *  —  >- 

•  x  — 


<0— W-— x— a_a 
o  —  —  x  u.  i  <m_  — 
a  m *  ♦  ♦  —  »—  a  — ■■  fsi 
O'  i—  —» —  y  3  •  a  a 
a  u.  —  X  O  •  -J  w*  a 
—  w-»Su  J  - - i  -J 


<N»  a  ♦  —a 
a  - — ■' —  r’l  X 

<\i  ■ 


Mi  yui  _ 

O  a  ~—  a 

o  O  a  a  »  «n  —  —  • 
S<  CT*  — XU.C 
*  — •  a  ■ — ■  a  a  <  m  a 
o.a'Oar—  ms 


— >wi  <v<  u.  m.  ac  — *  ♦ 

•ar  q  — -  w  m.  ■ — '  X  ■—  1 

X  a  »  o  3  — —  a  m'-’ 

♦  •— *  &  a  »  ct3  X  a 

—  00  a  —  a  -in  3  - 

<M  ^  X  mJ  •  a 

a  u»  3  •—  O  >—  —*  • 
a  M-S'-'V/la  JNa 

r\i  O’ 

■—»  O  V/7  M-  «  !M  J  —-MM 

a  jsowOa  ■nu.'-' 

—  X  «  O  •  a  —  u-  a 

^  ♦  3  vj  ’ —  O  1  ► 

b-a  a  x.*0'~"~ 'OX' 
m.  a  —  ♦  ~  a_  at :  ■ 

Mi  a  — ■  rsi  ^  X  *  CL. 

o  —  — ■  a  M.  «-  a.  a  • 

j  ">  t  •  Mi  a  ac 

X  b->M.  O  O  M.  ur»  O  a 

♦  m.  «j  >— » Q  m.  O  a  *>i 

IN4M.O"  XmJ^— >■ 
a  ^o  —  ♦  O  O  —  — 

a  O  X  cn  —  W  •  "N  ♦ 

a.  o  ♦  —  a  S  -—  —  O' 

a  x  mi  — — «  a  ac  u.  a 

■ — ■ — -  j  ib  kM/i  <  a  a  i' 

!\ja  mi  m.  m  _j  ►—  uj  •  • 

^TC’a  2  ^  U  x  a. 

*n-'U«  •  *—  X  “•  • 

Mi  —  O’  x  o  a 

5  — a  cn  •  x  j_ 
j  -—vl  •  — a  Q  — — 

X  *  X  —  i1«  U-  '— 

— •• — ■  <  N"  33  I.  « 

«  —a  a 


O’ 

a 

a 

a 

O 


Mi  —  —  a  —  «  « 
—  »"N  W  <»  x 


<  — 

—  Oa  a  —a  —  — — 
ji  o  -J  —  —  aa  Mi  •— 

—  Mt  a  'b.^Xa  ■ — ■  O-i 
Mi.x— — O  —  X«  a 

—  :m  ^  ji  j  n  a 


<j 

X 

3* 


lJ 


•m)  O  <-  ’ 


I 

J 


MUT-t 


SUBROUTINE  RUUOEK 
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IMW.L4.GT.O.)  GO  TO  1 

ENSURE.  TIHE  LAG  HAS  PASSED  BEFORE  RUOOER  MOVES. 


U(  T6.LE.TLAGM)  GO  TO  3 
If (UAbS(UlL4 ) . LT.DCLW)  W)  TO  J» 


If (fl 1.LI.TIAGR)  GO  TO  8 


So  J 

►—  i—  ■  o 

•  «  CM  —  **»  Si 


a 


SUBROUTINE  DEPTH 


y  tf(DLLB.GE.ttOWIAX)  DlLb-BO**AX 
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IF(UABS(UELB) .6T.ADEL1 )  GO  TO  11 
LHSUHL  MAX  PLANE  ANGLE  IS  NOT  EXCEEDED. 

B  IP(UABS(UELB).GE.bOWrtAX)  DELB-BQMHAX 

IF  WU  Ul VE/ASCENT  ANGLE  IS  EXCEEDED  UNITE  DIAGNOSTIC. 


IMATHtlA.Gi.MAJtAWG)  GO  TO  G 
U)  TO  91 


L 


O 

u. 

3 

3 


V* 

3 

3. 


* 


3 

_ i 
* 


3^ 

42 

33 


3 

»n 


o 

3 


co 

« 

5 


CO 

3 


X 

3 


3  I 
5  2 


•  -n 

o  — 

«  ^ 


< 

H- 

-J 

—  3 


9  -1 
.  3 


«  « 
3  3 


33 

*—  H- 
-4  — i 

33 


■ 

* 


<  « 


3* 

33 


55 


O 

3 

o  • 
o  « 

o  « 
x  • 


•  9 

«  (Vi 

«  * 


at 

:* 

o 


a  — 


X  - 
>-  O  • 
—  CM  • 
n»  * « 
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x*  CM  CM  .— 

*  x  5  o  5 

SSoooo^-o 

3  3  ►— 

►—  »—  o  Q  UJ 

uiuiQQXQo  ae 

ac  ac  <3  <3  o 


ooxoxo  •  ■ 


*-  < 
z  ►- 

o 

o  a. 

3 
<  X 

M  — 

O"  O'  O  (O  »--  O  O  K- 
ujuJOZOZOO 

aaaaaaaa 

oooooooo 

3333 3333 

oj  Gk. 

**•••••♦ 

00000330 

3  ^ 

.  55555555 

s 


oooooooo 
►—  ^ 

OoOojOoO 


_  _  —  — i  ~  —  —  —  !M 

a a  335535 3 


oooooooo^. 

355535553 


000003090 


a  a  cm 
a 


v-1 

< 


332 

--3 


3 

< 


X  3 
Z  u. 


3 

3 


3 

OR 


3OOOJOOOO 


aooaooooo 

MNNVMNNNW 


CM 


3  0  0 


o  o  o 


0(0  0 


o  o  o 


*  I 


uttu-^ 


UtLb— BOMMAX 


SUBROUTINE  STERN 


|F(DtLS.Ut.STL»U)  ULLS-STEWU 


lF(OELTb.LE. DELTA)  0ELS-A0CL3 
If (UllTB.GT. DELTA)  DUS-OELS-OEITA 
Go  TO  W 

tHSURL  T IMF  LAG  HAS  PASSED  BEFORE  PLANE  MOVES. 


ti  ni  oo  (swu*vr«i)n 

0-/| 


nnjKii-M  <m 

f»6  01  00 
0-STK1  SI 
66  01  00 

vm<HSTHi*smo  m 


nin*oii-oii  zi 
m  01  no 


■5 

4 


-J 
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8.3  IEQTTF 

A  -  Input  matrix  of  dimension  :i  by  N  containing  the  coefficient 

matrix  of  the  equation  AX  ■  B. 

On  output,  A  Is  replaced  by  the  LU  decomposition  of  a 
rowwise  permutation  of  A. 

H  -  number  of  right-hand  sides.  (Input) 

N  -  order  of  A  and  number  of  rows  In  B.  (Input) 

1A  -  number  of  rows  In  the  dimension  statement  for  A  and  B  In 

the  calling  program.  (Input) 

B  -  Input  matrix  of  dimension  M  by  *1  containing  right-hand 

sides  of  tne  equation  ax  •  B. 

On  output,  tr.e  H  by  '*■  solution  t  replaces  3. 

IDGT  -  Input  option. 

If  IDGT  Is  greater  than  J  the  elements  of  A  and  B  are 
assumed  to  be  correct  to  I06T  Decimal  digits  and  tne 
routine  performs  an  accuracy  test. 

If  IDGT  equals  zero,  tne  accuracy  test  Is  bypassed. 

■KARLA  -  wort  area  of  dimension  greater  than  or  equal  to  *1. 

ItR  -  error  parameter 

terminal  error  ■  123  ♦  N 

H  •  i  indicates  tn«i  A  is  algori tnsical !y  singular. 

warning  error  ■  32  ♦  1 . 

N  •  2  Indicates  that  tne  accuracy  test  failed.  The 
computed  solution  may  be  In  error  by  more  than  can  be 
accounted  '"or  ay  tie  uncertainty  of  trie  data. 

CALL  IE3T1F  (A,  A,  IA,  B,  IDGT,  WKAREA,  ICR) 


